Telegraph-type versus diffusion-type models of turbulent relative dispersion 
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Properties of two equations describing the evolution of the probability density function (PDF) of 
the relative dispersion in turbulent flow are compared by investigating their solutions: the Richard- 
son diffusion equation with the drift term and the self-similar telegraph equation derived by Oga- 
sawara and Toh [J. Phys. Soc. Jpn. 75, 083401 (2006)]. The solution of the self-similar telegraph 
equation vanishes at a finite point, which represents persistent separation of a particle pair, while 
that of the Richardson equation extends infinitely just after the initial time. Each equation has a 
similarity solution, which is found to be an asymptotic solution of the initial value problem. The 
QQ ' time lag has a dominant effect on the relaxation process into the similarity solution. The approach- 

' ing time to the similarity solution can be reduced by advancing the time of the similarity solution 

, appropriately. Batchelor scaling, a scaling law relevant to initial separation, is observed only for the 

■ telegraph case. For both models, we estimate the Richardson constant, based on their similarity 
» I ' solutions. 
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' The relative dispersion of particle pairs is fundamental to the turbulent research and has many practical applications 
1—1, to environmental and industrial problems, such as the transport of pollutants in the atmosphere and the fuel mixing in 
' engines. The study of turbulent relative dispersion has a long history since the pioneering work by Richardson^, who 
, observed the anomalous dispersion or superdiffusion of particle pairs and proposed the diffusion equation describing 

• ■ the evolution of the probability density function (PDF) for the pair separation. Several attempts to modify the 
', Richardson model have been presented^ii^i^. Recent experiments and direct numerical simulations (DNSs) showed 

■ that the Richardson model can well represent the separation PDF in the inertial subrange^s^ii'^'^. Nevertheless, the 
best model of the relative dispersion has not yet been determined decisively. This is because the inertial ranges 
achieved both in experiments and DNSs up to now are not wide enough to observe the full superdiffusive behavior. 

CN • The relative dispersion in turbulent flows has often been modeled with the Reynolds-number dependence and the 
contribution from the three ranges of turbulence, namely the energy-dissipative, the inertial and the energy-containing 
rangesiSiiiii^. Indeed, the unified description of the relative dispersion is important especially for comparison of data 
among models, experiments, and DNSs at moderate Reynolds number. However, we limit our attention to the 
relative dispersion in the infinitely-extended inertial range, namely in the case that Re ^ cxj, where the superdiffusive 

• ' behavior can be observed, and discuss the governing equation of the separation PDF there. As is already mentioned, 
I , the Richardson model is supported by several experiments and DNSs. However, there are a few shortcomings in this 

J — ' model. One of those of the Richardson model is the absence of the long-time correlation of the Lagrangian relative 
, velocity of a particle pair— li^, which exists in real turbulent flows. This correlation gives rise to persistent separation 
of pairs, which should obey the self-similarity as shown in the recent DNS^^. 

In order to overcome this shortcoming, Ogasawara and Tobi^ devised another model of turbulent relative dispersion, 
the self-similar telegraph model. This model is derived taking into account the persistency of pair separation on the 
' basis of Sokolov's pictur o^^d^ , where the relative separation consists of persistent expansion and compression. As a 
result, the governing equation of the PDF of relative separation in this model has a second-order time derivative term. 
This is the reason why this model is called 'telegraph,' while in the other previous theories such as the Richardson 
model the governing equation of PDF is diffusion- type. The term 'self-similar' comes from scale-dependent coefficients 
reflecting the self-similarity of the dispersion process in the inertial range. Owing to this inclusion of the self-similarity, 
solutions of the self-similar telegraph equation do not approach those of its corresponding diffusion equation, whereas 
those of the usual telegraph equation appearing in the problems of molecular diffusion or Brownian motion quickly 
relax into those of the diffusion equation. This indicates that the self-similar telegraph model has essentially distinct 
properties from those of its diffusion- type counterpart. 

In this paper, we investigate time-integrated solutions of the self-similar telegraph equation as well as those of its 
diffusion-type counterpart. In the literature, the realizability of the similarity solution of the Richardson diffusion 
equation has implicitly been assumed, and thereby the agreements between the similarity solution and the separation 
PDF obtained by experiments or DNSs have been discussed. However, the similarity solution becomes a delta function 
at the origin in the limit t — > 0. This requires that a particle pair be located at the same place aXt — Q. This condition is 
difficult to set in experiments and DNSs. In addition, the viscous effects in the dissipation range may also contaminate 
the pure inertial-range behavior of the dispersion process. For this reason, the initial separation of a pair is finite in 
experiments and DNS, and hence the initial condition with a nonzero separation is desired. 
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Furthermore, the recent experiment^ - suggested that the separation PDF obtained by their experiment agreed well 
with the similarity solution of Batchelor's diffusion equation rather than Richardson's mentioned above. However, their 
claim is contradictory in that the similarity solution of the Batchelor diffusion equation implies the usual Richardson 
scaling law for the mean-square separation (r^) cx t^, whereas they found the (modified) Batchelor scaling law instead 
of it, as was also reported in Ref. 18. These results stimulate us to investigate the short time behavior of the solution 
of the Richardson diffusion equation with finite initial separation. To our knowledge, this issue has never been treated 
in the literature. 

Motivated by the above, we numerically solve the governing equations of the PDF under appropriate initial condi- 
tions with finite separations, and observe the behaviors of the time-integrated solutions. Comparisons of the solutions 
between the telegraph-type and the diffusion-type model are also made to reveal the characteristics of the models. 

The remainder of the paper is organized as follows. In Sec. [IT] we introduce the self-similar telegraph equation 
together with its diffusion- type counterpart. Section Hill explains the settings of the simulation. The results of the 
simulation is presented in Sec. lIVl where the Richardson constants are also obtained and compared for the two models. 
In Appendix El the decrement of the total probability is examined for the self-similar telegraph model, because the 
conservation of it is not guaranteed for this model. 



II. SELF-SIMILAR TELEGRAPH AND PALM EQUATIONS 



In this section, we introduce the self-similar telegraph model of turbulent relative dispersion. In this model, the 
evolution of the spherically symmetric PDF of pair separation P(r, t) in d-dimensional isotropic turbulence is described 
by the following self-similar telegraph equation^^ 
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where Tc(r), D{r) and v{r) are a characteristic time scale, a diffusion coefficient and the Lagrangian relative velocity, 
respectively, of a particle pair of separation r. The parameters A and a characterize the turbulent field considered. 
The coefficients Tc{r), D{r) and v{r) are assumed to obey the following scaling laws: 
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Here, A is a dimensional constant and s a scaling exponent: s — 2/3 for Kolmogorov scaling and s = 2/5 for Bolgiano- 
Obukhov scaling. The last term on the r.h.s. of Eq. ([T]) is a drift term, with the drift velocity being —(jv{r). A~^ 
represents the persistency of the separation, and corresponds to the persistent parameter introduced by Sokolov^^. 

In the derivation of Eq. ([T])^, a parameter 5 was also used, and has a relation with A and a via \a = d — 2s + 6 . 
The physical meaning of 5 is the difference between the two transition rates of the direction of the separation, from 
expansion to compression and from compression to expansion. Hereafter, we use (5 as a control parameter instead of 
cr, following the previous paper a^^'^^ . 

If the effects of the finite separation and the finite correlation of the relative velocity are not considered, the first 
term on the l.h.s. of Eq. ([!]) is omitted, and then Eq. ^ can be reduced to the following diffusion equation: 
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This equation has the same form as the Richardson diffusion equation with the drift term. If cr = 0, i.e. S = 2s — d, 
Eq. is the Richardson diffusion equation itself. The addition of the drift term to the Richardson equation was 
discussed by Palm in 1957^*^. Recently, Goto and Vassilicos^^ derived the same equation as Eq. We refer to 
Eq. ^ as Palm equation. 

Before we perform the time integration of Eqs. ([1} and ([3]), the nondimensionalization is made for convenience. We 
normalize space by the nonzero initial relative separation R 0) of particle pairs and time by the corresponding 
timescale Tc{R). Consequently, new dimensionless space and time variables are defined as r ~ r/R and t = t/Tc{R) = 
At/R'^ respectively. We also introduce the normalized separation PDF P{r,i) = RP{r,t). Substituting Eqs. ^ into 
Eqs. IJ) and Q, eliminating the dimensional quantities and dropping the superposed tildes from the dimensionless 
quantities give the following nondimensionalized versions of the equations: 
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In the following, we treat Eqs. (|4]) and (O instead of Eqs. ([T]) and ([3]). 



III. SIMULATIONS 
A. Parameters 

We set the values of the parameters d and s to 2 and 2/5, respectively, assuming the case of two-dimensional free 
convection (2DFC) turbulence, which is characterized by the Bolgiano-Obukhov scalin g^°i^^i^'^ . This is because we 
have recently estimated the values of the control parameters, A and S, using exit time statistics^^. For A and S we 
took the values 5.2 and —0.77, respectively, from Ref . flQl. Since one of the main goals is to compare the properties of 
the different models of turbulent relative dispersion, this choice of the parameters is sufficient unless the variation of 
the parameters yields qualitative changes of the solutions of the model equations. In addition, it is worth noting that 
the one-particle dispersion in the 2DFC turbulence was also investigated only recently^^ 

B. Initial and boundary conditions 

We aim to examine the case that the pairs of the particles have the same nonzero separation at the initial time, 
and therefore, the initial condition should be represented by a Dirac's delta function located at a certain finite point. 
To mimic this initial condition, we employed the following finite-width PDF as an alternative initial condition: 

P(r,0)=|i(l + ^°^('^^)) l^-l|<" (6) 
[ |r-l|>w, 

where w represents the width of the initial PDF. As is mentioned in the last part of Sec. [TTl the spatial scale is 
normalized by the initial separation, so that the initial PDF ([6]) is distributed around 1. Hereafter, we regard 1 as 
the typical initial relative separation. Note that this value cannot be strictly referred to as the mean value since the 
PDF should be spherically integrated to obtain the mean value of relative separation in two or three dimension. In 
the following, we consider the two cases, w = 1 and 0.1, to see the effects of the width of the initial PDF. Narrower 
widths than 0.1 cause the numerical instability of the scheme for the present interval of grids and steeper ridges at 
the bounds of the time-integrated solution of the self-similar telegraph equation (see Subsec. HV Al and Fig. [1]). 

For the self-similar telegraph equation another initial condition [dP{r,t)/dt]t=o = is required. This condition 
implies the symmetry between the extending and compressing pairs at the initial time, which reflects the situation 
where the initial placement of pairs is uncorrelated with the velocity field of turbulent flow. Finally, we set the 
boundary conditions to P(0,t) = P{oo,t) = 0. 

C. Numerical method 

Time integration of Eqs. ^ and (O is performed with the Crank-Nicolson scheme, which is second order accurate 
in time, after the transformation of the spatial variable. We define a new spatial coordinate r' as r" /s and r/ {r + L) 
for Eqs. (HI) and (O, respectively. Here L is an arbitrary parameter, which determines "local density" of grids in r 
space. The former transformation leads to the constant Courant number c = r^~*(Af/Ar) = At/Ar' with uniform 
grids in r' space, while the latter transforms r € [0, oo) to the finite computational plane r' G [0, 1]. 

For the self-similar telegraph model, grids are added above the maximum grid of the computational plane and 
remove the grids of the same number as the added ones when the maximum separation of the PDF nearly reaches 
the maximum grid. In this way, the number of the grids is always conserved during the computation (although the 
interval between the grids are increased at every change in the grids). The grids of 82, 000 and 100, 000 were used to 
compute the solutions of Eqs. (jl]) and respectively. For the self-similar telegraph model, we set the initial grid 
interval Ar^^^^ to 10~*, while for the Palm model, we set L to 0.2 or 10 for the short-time and 100 or 1000 for the 
long-time behavior. 
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FIG. 1: Time-integrated solutions of the self-similar telegraph equation Q and of the Palm equation ((5}. The solid lines denote 
the case of w = 1, while the dashed lines w = 0.1. {a.) t — 1 for the telegraph case, (b) t — 1 for the Palm case and (c) t — 4 
for the telegraph case. The insets are the same plots in linear scale. 



IV. RESULTS AND DISCUSSION 



A. Short-time behavior 



The time-integrated solutions of Eqs. ^ and ([5]) are shown in Fig.[T]at different times. While the solution of Eq. ([5]) 
extends infinitely at finite time, that of Eq. ^ bounds at a finite point. This point corresponds to the maximum 
separation of pairs, which is exactly obtained by direct integration of the equation dr/dt — v{r) with the scaling law 
of the velocity difference (Pc)) to give in the nondimensionalized form 

w(^) = (C(0) + s^)l/^ (7) 

The existence of the maximum point in the solution is a manifestation of the finite separation of a pair, which is 
included in the self-similar telegraph model. 

For the case oi w = 0.1 (Figs. [BJa) and[IJc)) of the self-similar telegraph model, we can observe ridges at the 
edges of the PDF. The similar behavior has been found in the model of Ref. |U based on a Levy-walk stochastic 
approach. The PDF predicted by their model has peaks at the side edges. Since our approach is similar to theirs, 
the reason for the presence of the ridges should be consistent: the distribution of the relative velocity is not taken 
into account. Particle pairs having never been compressed or expanded from the initial time should be accumulated 
at these advancing or receding ridges, respectively. To see this, let rmaxiO) in Eq. ([7]) be equal to 1, the center of the 
initial PDF. Then, r„iax{'^) — 2.32, consistent with the position of the advancing ridge in Fig. [TJa). However, these 
ridges get smaller as time elapses owing to the effect of diffusion, i.e. the random changes of the direction of relative 
velocity. In Fig. mc), at i = 4, the receding ridge disappears and only the remnant of the reduced advancing ridge 
can be observed in the case of ui ~ 0.1. 
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FIG. 2: Similarity solutions of the self-similar telegraph equation (j4| (solid line) and the Palm equation ([5| (dashed line) in 
the similarity form F(r]) = {t/X)^^" Ps{r,i). The similarity variable rj is normalized by r]c = s\, corresponding to the maximum 
separation for the self-similar telegraph equation Q. s = 2/5, d = 2, \ = 5.2 and 5 = —0.77. The inset is a linear plot of the 
same figure. 



B. Long-time behavior 

Equation Q can be reduced to a second-order ordinary differential equation with the similarity variable rj = Xr'^ /t. 
Therefore it has a similarity solution, and so for Eq. (O. Although the similarity solution of Eq. ^ has no analytical 
form, that of Eq. ([5]) does and reads 

''.(M) = ^($;)" »P (-$;). (8) 

where C is the dimensionless normalization factor 

and T{x) = e-H^'-'^dt is the gamma function. The subscript s of Pg denotes the similarity solution. In Fig. [2l we 
show the similarity solutions of Eqs. (|4|) and ([5]) in the similarity form F{ri) = [t/ X)^/^ Ps{r, t), with the values of the 
parameters specified in Subsec. IIII Al 

In contrast to the solution of the Palm equation, that of the self-similar telegraph equation has an upper bound 
in space, and vanishes for rj > rjc = sX, which is attributed to the singular point of the reduced ordinary differential 
equation^^. This yields the maximum separation Vmaxit) = (si)^^^, which is also obtained by neglecting rmax{0) in 
Eq. 0. 

The similarity solutions tend to a delta function at the origin in the limit t —>■ 0^ i.e. Ps(^, 0) = S{r), as is expected 
by Eq. ([8]) and the fact that rmax{0) = for the similarity solution of the self-similar telegraph equation. This initial 
condition is different from that adopted in the present paper. However, we inferred that the effect of the initial 
separation on the PDF is negligible after a long time, and thus the solutions of Eqs. ^ and (O starting from our 
initial condition ([5]) approach the corresponding similarity solutions. In order to confirm this, we further advanced 
the time-integrated solutions. The resulting PDFs are transformed into the similarity form F{ri) and are plotted in 
Fig. [21 From this figure, we found that the solutions with the initial condition ^ approach the similarity ones as 
time elapses, and finally almost reach them at t = 4092, so that we cannot distinguish them. This suggests that the 
similarity solutions are asymptotic solutions of the initial value problem. The same behaviors as that shown in Fig. [3] 
can be obtained for w = 0.1. 

In order to analyze this behavior quantitatively, we introduce the root-squared difference between the time- 
integrated and the similarity PDF as 



\\SP\\ = JI \Pir,t)-P,ir,tfdr. (10) 



This quantity vanishes if the separation PDF has the same form as that of the similarity solution. Thus, we can 
quantify the degree of the proximity of the two solutions by calculating this quantity. 
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FIG. 3: Time-integrated solutions of (a) the self-similar telegraph equation (|4]) and (b) the Palm equation ((5]) in the similarity 
form F{r]) = (t/Xy^" P{r,t), under the initial condition ^ with uj = 1, at time t = 12, 60, 252, 1020 and 4092 (from bottom 
to top at the maximum point), along with each similarity solution (solid line). The insets are linear plots of the same figures. 




FIG. 4: Temporal evolution of the root-squared differences between the time-integrated and the similarity solutions, \\5P\\ 
or ||(5P||t;, defined in (|10|l or (|15|) . for both equations and widths. The upper curves denote while the lower ||(5P||ti: 

advancing times tat in Eq. (|15p are set to 8.5, 8.3, 8.3 and 8.2 for the telegraph with w = 0.1, the telegraph with it; = 1, the 
Palm with w = 0.1 and the Palm with «; = 1, respectively. 



In Fig. m we plotted the temporal evolution of \\6P\\ for both models and widths. From this figure, the scaling 
laws \\SP\\ cx can be found in both short-time and long-time regimes with (3 ~ 1.1 and 2.2, respectively. The 
transition between the two regimes occurs around t = 10. No qualitative difference between the two models is found. 
The difference between the two values of w in the initial condition ([6|) can scarcely be seen throughout the whole 
time scales shown in the figures. In this case, the relaxation processes into the similarity solutions appear to be 
rather universal, independent of both models and widths. However, as is illustrated below, these apparent universal 
behaviors are mainly due to the effect of time lag. 

As an explanation of the origin of the above scaling behaviors, we consider the effect of time lag, because the 
similarity solutions are localized at the origin in the limit t — > 0, whereas the initial condition for the time-integrated 
solutions ([6]) implies nonzero initial separations. We assume that the initial condition ([6]) can be approximately 
described by the similarity solution at a certain time to- Then, its subsequent time evolution is the same as the 
original similarity solution after to only for the Palm case, since for the telegraph case it depends on the additional 
condition [dP{r, t)/dt]t=o. We can now replace P(r, t) by Ps{r, t + to) in the right hand side of Eq. (fTU]) . and substitute 
Eq. dl]) to obtain 

\Ps{r,t + to)-Ps{r, i)|'rfr 
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FIG. 5: Fitting of \\SP\\ with the root of Eq. ([TTjl . in the case of the Palm model ^ with w = 1. The solid and the dashed 
lines represent \\5P\\ and the root of Eq. with the fitting parameter to — 8.2, respectively. The dotted and the dot-dashed 
lines are the roots of Eqs. (|13|l and (|14|) . respectively, with to — 8.2. 



where 

4s - 2(5 - 1 



(12) 



We examine the two limiting cases for Eq. (jlip . namely the long-time and the short-time limits. The calculation for 
the former case is straightforward and gives 



\p,ir,t + to)-Ps{r,ttdr^-^(Q ' r(7)(7S^ + 1) f^V (^o « 0- (13) 







For the latter case, it is obvious that 



(J2 / ^2 \ i/s 



\Ps{r,t + t,)^Ps{r,t)\^dr^ P^{r,t)dr ^ j:^^ i^- j r{j) (io » 0, (14) 

since only the P^{r,t) term in the integrand of Eq. pT|) diverges in the limit t — > 0. Equations and (fH|) lead to 
the scaling laws = t~^'^^ and t^i/(2s) — t^^-^s^ respectively, for the root-squared difference. This explains 

the above scaling behaviors of \\SP\\ in FigHl 

Figure [5] displays the result of the fitting of \\6P\\ in the case of the Palm model ^ with w = 1, shown in Fig. IH 
with the square root of Eq. pT|) . Here, the time lag to of value about 8.2 for the square root of Eq. pT|) gives the 
best fit of The agreement is excellent. This result ensures our expectation that the time lag has the dominant 

effect on the decay of \\6P\\ in time. 

By removing the effect of the time lag, more rapid approach of the time-integrated solution to the similarity solution 
can be realized. In Fig. [5J we also plotted the following quantity: 



|<5P||t, = ^y^ \Psir,t + tat)-P{r,t)\''dr. (15) 

Here, tat is an advancing time of the similarity solution, which is chosen such that ||(5P||t; has the minimum value 
around t ~ 1000. Specifically, tat — 8.5, 8.3, 8.3 and 8.2 for the telegraph model with w = 0.1, the telegraph model 
with w — 1, the Palm model with w — 0.1 and the Palm model with w = 1, respectively. These values are consistent 
with the above optimal value of to in Eq. pT|) . resulting from the fitting (Fig.O- As is obvious from the figure, 
is reduced by two or three orders of magnitude compared with ||(5P||. Moreover, the algebraic decrease is no longer 
seen for both short and long times. 

In Fig. \6\ we show the approach of the time-integrated solutions to the corresponding similarity solutions advanced 
by tat — 8.3 and 8.2 for the self-similar telegraph and the Palm model, respectively. From the comparison between 
Figs. [3] and [6l it is obvious that the time-integrated solutions approach the similarity solutions with the time lag much 
faster than those without the time lag, as is expected from Fig. ID At t = 252 already, the time-integrated solutions 
almost accord with the similarity solutions. Hence, the approaching time to the similarity solution becomes shorter 
than that without the time lag as in Fig. [31 by one order of magnitude. 
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FIG. 6: Same as Fig. |3]but for another variable r\ = (Ar'')/(t + tat) with tat ~ 8.3 and 8.2 for the self-similar telegraph and 
the Palm model, respectively. Here, the solutions a.t t = 1020 and 4092 are omitted, since they are almost indistinguishable 
from the corresponding similarity solution. 
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FIG. 7: Short-time behaviors of the mean-square separations. Here, their initial values, {r^(0)) = 1 -f (1/3 — 2/it^)w^, are 
subtracted. The solid lines are for w — 1 and the dashed ones for w = 0.1. The straight dot-dashed lines represent 1.3t and 
2.4:t^ for the Palm and the self-similar telegraph model, respectively, for reference. 

C. Batchelor scaling 

For times shorter than the characteristic time of an eddy of size of the initial separation, the difference between 
the mean-square separation and its initial value is expected to grow in time as t^, reflecting the initial ballistic or 
persistent motion of particle. This scaling law was first derived by Batchelor—. and we refer to it as Batchelor scaling. 

To see this Batchelor scaling for the present time-integrated solutions, we plot in Fig. [7] {r'^{t)) — (r^(0)) versus 
time for the above solutions. For the self-similar telegraph model this scaling law is clear. This result implies that 
the time derivative of the mean-square separation vanishes at the initial time for the self-similar telegraph model, 
which reminds us of the initial condition [dP{r,t)/dt]t=o = 0. However, the scaling law proportional to t is observed 
for the Palm model. This arises from its incapability to represent the condition [d{r'^{t))/dt]t=o = 0, or equivalently 
[dP{r,t)/dt]t=o = 0. Therefore, the Palm model cannot satisfy the initial symmetry condition of pair separation 
mentioned in Subsec. IIII Bl Notice that from Fig. [7] the behaviors for the two different initial widths are almost the 
same, and hence it is concluded that the presence of the ridges only in the case of w = 0.1 for the telegraph model at 
short times, shown in Figs.[Ua) and[ljc), does not affect this scaling law. 

In the previous subsection, we showed that there exists the optimal time of the similarity solution for minimizing 
the relaxation time. This may lead to another scaling law of the mean-square separation. If we assume that the 
root mean-square separation obeys the same temporal evolution as that of the maximum separation ([7]) except the 
coefficient of the time, we may write the mean-square separation as 

(r2(i)) = ((r2(0)>^/2+rf)2A, (16) 
where c is a constant. This expression was also proposed by Goto and Vassilicos^-'^ using the concept of the doubling 
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FIG. 8: Comparison of the scaling law (|16p with the results of the simulations. Panel (a) corresponds to ct and (b) to c in 
Eq. ([16}. The straight dashed line in (a) represents Q.2t for reference, w = 1 for (a). The solid and the dashed line in (b) 
denote the cases of to = 1 and 0.1, respectively, for the Palm model. 



time, and was adopted in the recent paper—. Equation (fT6|) can be interpreted as the Richardson law with the 
advancing time, {r^(t)) cx (tat + t)'^^'^, if we regard (r-^(O))'''^^ as ctat- This means that the scaling law (fT6|) exactly 
holds for the similarity solution advanced by tgt In fact, Eq. (fT6|) was employed to determine the Richardson constant 
in DNS with finite initial separations in Ref. 21. 

In Fig. [S^a) , this scaling law is observed for the Palm model. Note that this scaling law tells nothing about the 
Batchelor scaling. In fact, as is seen in Fig. [5]Ja), Eq. is not satisfied by the self-similar telegraph model, which 
satisfies the Batchelor scaling. 

Figure[8jb) shows the time dependence of the coefficient c in the right hand side of Eq. (I16p for the Palm model with 
w = 1 and 0.1. The values of c slightly vary in the region 1 < i < 10, which means that Eq. (fT6|) cannot exactly capture 
the behavior of the mean-square separations of the time-integrated solutions in the whole time regime. Moreover, the 
basic assumption (r^(0))^/^ ~ ctat and the estimated value of the advancing time tat in the previous subsection lead 
to the value of about 0.12 for c for each width of the initial distribution of pair separation. This value conflicts with 
that obtained in Fig. [8jb). Therefore, the advancing time tat of the similarity solution, introduced in the previous 
subsection, does not have a significant meaning for the short time behavior. The variation of c is smaller for w — 1, 
implying that Eq. (|16[) is more appropriate for the smoother initial PDF, and becomes less appropriate as the initial 
PDF gets sharper and closer to the shape of the delta function. In the long-time limit this value is related to the 
Richardson constant, which is the subject of the following subsections. 



D. Richardson scaling 

For times much longer than the timescale of the initial separation and much shorter than the integral timescale, 
the well-known Richardson t^ law of the mean-square separation is expected to hold in SDNS or 2DIC turbulence 
governed by the Kolmogorov scaling. In the Bolgiano-Obukhov scaling, however, the mean-square separation should 
grow as t^ . This scaling law is clearly seen for the two models in Fig. [5] for the long time. 

However, now that we previously showed the agreement of the time- integrated solutions with the similarity solutions 
in the long-time limit, we can use the similarity solutions to see that scaling law. From the similarity form of the 
PDF F{ri) = (t/Xy/Tsir^t), it is easy to derive 

{r^)=(^iy^l jn^l^-^F{^)dr^. (17) 

Here, the integration on the right hand side is regarded as a constant for the similarity solutions, so that we can 
rewrite Eq. pT|) as 



2/s 



r) = G[j] ■ (18) 



Thus we can find the t^ dependence of the mean-square separation in the Bolgiano-Obukhov scaling s = 2/5. 
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FIG. 10: The coefBcient of the Richardson i}^" law (|18|) for the similarity solutions, G, as a function of the parameters A and 
(5. The upper dashed lines refer to the Palm model, Eq. (|19p . and the lower solid lines the self-similar telegraph model, (a): 
s = 2/5 for the Bolgiano-Obukhov scaling, (b): s = 2/3 for the Kolmogorov scaling. 



E. Richardson constant 



One of the main interests in turbulent relative dispersion is the determination of the universal coefRcient of the 
Richardson law, namely the Richardson constant''. From Eq. (US]), G{C ^ / X)^/^ is regarded as the Richardson constant, 
where is the nondimensional part of the dimensional constant A (see Appendix [B|) . 

First, we determine the value of G. To do this, we have only to calculate the integration in Eq. p7|) for the similarity 
solution. For the Palm model, this can be easily carried out using the similarity solution ([5]), and in terms of the 
gamma function we write it as 

Note that this value is independent of the value of A. This is because the Palm equation ([5]) becomes independent of 
A if we use rescaled time t — t/X. 

For the self-similar telegraph case, however, the dependence of A cannot be removed by this rescaled time, and hence 
G depends on A. Although we do not have the explicit expression of G for this case, we can evaluate it numerically, 
and the evaluated value are shown in Fig. [TO] for various realistic values of A and S with those of the Palm model. 
For s = 2/5, A = 5.2 and 6 — —0.77, we have G = 0.34 and 0.66 for the self-similar telegraph and the Palm model, 
respectively. 

The value of G for the self-similar telegraph case accords with that for the Palm case in the limit of A ^ oo, since 
Eq. (HI with the rescaled time i tends to the Palm equation for this limit. As can be seen from Fig. [TO] G of the Palm 
model is always larger than that of the self-similar telegraph. This means that the inclusion of the effect of persistent 
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TABLE I: Estimated values of the Richardson constant and the relevant parameters. The dashes denote unavailable values 
from the papers. Here, the values are estimated by us, except for the rightmost column. For the estimation of A and S, see 
Ref. [H. 







parameters 




Richardson constant from 




paper 


s 




A 


S 


the telegraph 


the Palm 


the original 












model ID) 


model (0 


paper 


Boffetta and Sokolo^*^^ 


2/3 


0.72 


14 


-1.48 


1.2 


4.4 


3.8 


Goto and Vassilicos^i 


2/3 


0.84 




-0.87 




3.9 


6.9 


Ogasawara and Toh^^ 


2/5 


0.62 


5.2 


-0.77 


0.030 


0.058 





separation for the model suppresses the relative dispersion. 

Now, wc calculate the Richardson constant for the two models. The remaining component of the Richardson 
constant, C^/A, must be determined from the characteristics of turbulent flows. This factor is estimated in the way 
summarized in the Appendix [Bj Note that the value of A is also required for the self-similar telegraph model, since G 
depends on A in this case. For our previous DNS of 2DFC turbulence, we have estimated the Richardson constant as 
0.030 and 0.058 for the self-similar telegraph and the Palm model, respectively. These values are smaller than those 
estimated of the Richardson constant for the 2DIC case, also shown in Table [H by two orders of magnitude. This 
might be one of the peculiarities of relative dispersion in 2DFC turbulence. From the upper two cases in Table [H 
the Richardson constants estimated using the Palm model appear closer to the values obtained by their DNSs (the 
rightmost column). However, in their original papers^*"^^, additional assumptions were imposed to obtain their values, 
shown in the rightmost column of Table [J the satisfaction of the Richardson diffusion equation for Ref. and the 
scaling law of Eq. (|16p for Ref. [2l|, both of which do not hold for the telegraph model. Therefore, their values cannot 
be used to determine which model yields a more correct value of the Richardson constant. 

V. CONCLUDING REMARKS 

We have numerically solved the two different equations describing the temporal evolution of the PDF of the sep- 
aration of a particle pair in the inertial range of homogeneous and isotropic turbulence. In the simulations, 2DFC 
turbulence case is dealt with, because we have the values of the control parameters for that casein. The time-integrated 
solutions of both equations have been compared to characterize the models. There, the two initial conditions of the 
different widths of the initial PDF are imposed. However, the difference between the two widths only affects the 
short-time behavior of the PDF for the self-similar telegraph model. 

The self-similar telegraph model represents the finiteness of the separation and its persistency from the bounds of 
PDFs or the maximum relative separations, appeared in Figs. [T] and [H whose behavior is not seen in the diffusion-type 
counterpart (the Palm model) . The inclusion of the effect of persistent separation in the self-similar telegraph model 
would have an advantage in describing the initial persistent separation, which was recently found to strongly affect 
the relative dispersion for a relatively long time in experimentsiii^. 

In the case of the self-similar telegraph model, when the width of the initial condition ([6]) gets smaller, the ridges 
appear at the bounds of the PDF in the short-time regime. The uniqueness of the relative velocity is responsible 
for this, as is also shown in the model of Ref. [13 based on the similar approach to ours. However, in the real 
turbulent dispersion, the relative velocity has a distribution, and therefore these ridges cannot be observed in DNS 
and experiments. We will incorporate the distribution of the relative velocity into our model in future work. 

The long-time behaviors of the time-integrated solutions are characterized by a relaxation process into the cor- 
responding similarity solutions. It is found from the numerical simulations that the similarity solutions of both 
equations, ^ and ([5]), are asymptotic solutions of the time-integrated solutions. Thus, the realizability of the similar- 
ity solutions was corroborated, aside from the very long approaching time. We also investigated the decay rate of the 
root-squared difference between the time-integrated and the similarity solutions, II^J^II- The dominant contribution 
to \\SP\\ is the time lag between the two solutions. This contribution seems to be independent of both models and 
widths of the initial PDF. 

If we compare the time- integrated solution with the appropriately-advanced similarity solution, more rapid approach 
to it can be found. This is because the similarity solution is located at the origin at the initial time, while our initial 
condition ^ allows for nonzero initial separations. This result may give suggestions to the observation of the actual 
relaxation process of separation PDF in experiments or DNSs. 
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FIG. 11: Temporal evolution of the total probability, S, for various values of the parameters A and 6. (a): S — —0.77 and 
A = 3.5, 5.2, 6.8 and 8.3 from bottom to top. (b): A — 5.2 and S = —1.2, —0.77, —0.4 and from bottom to top. w = 1 for 
both. The values remain constant after t — 2. 



The Batchelor scaling law for the mean-square separation holds only for the self-similar telegraph model, while both 
models satisfy the Richardson scaling law. This difference for the short-time behavior between the two models is due 
to their capacities to satisfy the condition [d{r'^{t)) /dt]t=o — 0, or [dP{r,t)/dt]t=o — 0. 

On the other hand, the behavior of the mean-square separation for the Palm model can be well described by the 
Richardson scaling law with the time lag, Eq. ^TE\i . to some extent. However, as is suggested by the fact that the 
Batchelor scaling cannot be derived from Eq. (|16p. this scaling law neglects the effect of initial persistent separation 
of a particle pair. 

We have also estimated the Richardson constant for both models making use of the data of DNSs to determine 
the values of the parameters of the models. The Richardson constant for the self-similar telegraph model is generally 
smaller than that for the Palm model. Furthermore, both models predict smaller values for the 2DFC case than those 
predicted for the 2DIC case, by a factor of 10~^. This result should be confirmed in future DNS of 2DFC turbulence. 

Although we have dealt with the 2DFC turbulence, we have not yet observed the clear Richardson law in that 
case^^. We will attempt to achieve the Richardson law as well as the separation PDF in future work to determine 
which model is better for describing the relative dispersion. 



Acknowledgments 

This work was supported by the Grant-in-Aid for the 21st Century COE "Center for Diversity and Universality 
in Physics" from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan. Some of 
the numerical computations in this work were carried out on NEC SX-8 at the Yukawa Institute Computer Facility. 
K.K. thanks T. Matsumoto for suggestions on the numerical method and for invaluable discussions. S.T. was partly 
supported by the Grant-in-Aid for Scientific Research (C) (18540373) from Japan Society for the Promotion of Science. 



APPENDIX A: LOSS OF THE TOTAL PROBABILITY FOR SELF-SIMILAR TELEGRAPH MODEL 

Although for the Palm model the total probability S = j P{r,t)dr is conserved, the self-similar telegraph model 
does not assure the conservation of S. Thus we calculated S using the time-integrated solution and showed its 
temporal evolution in Fig. 1111 varying the values of the parameters A and d. We can easily see from this figure that S 
monotonically decreases for t < 2 and is conserved afterward. The variation of the value of w does not change these 
graphs, even quantitatively. Therefore, we may regard these results as universal, independent of the initial condition. 
Note that the loss of S is not the result of errors in accuracy of the numerical scheme used, since finer grids and time 
steps make the same figures. As is mentioned in our previous paper, the similarity solution fulfills the probability 
conservation. However, we cannot conclude from this that the time-integrated solution already has the similarity form 
for t > 2. 

We also plotted the dependence on the parameters A and 5 of the loss of the total probability, dS = 1 — S, a,t 
t = 2 m Fig. [121 There the scaling laws 0.44A~^'* and —0.013(5 -I- 0.011 are seen. Its decrease with increasing A is 
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FIG. 12: Dependence on the parameters A and 5 of the loss of the total probability, dS, at t = 2. (a): 5 — —0.77, (b): A = 5.2. 
The lines are fits to the data and represent 0.44A~^'* and —0.013(5 + 0.011 for (a) and (b), respectively, w = 1 for both. 



obvious from the fact that the self-similar telegraph equation ([4]) tends to the Palm equation fS]) as A ^ oo, keeping 
the rescaled time t appeared in Subsec. flV El unchanged. 

It is not difficult to make the self-similar telegraph model to conserve the total probability by reformulating it in a 
conserved form. The results will be reported in future papers. 



APPENDIX B: ESTIMATE OF C^/X 

The dimensional constant A can be decomposed into the nondimensional and dimensional parts as C^e^^^ and 

C^Sg^^ (ag)'^^^ for the Kolmogorov and Bolgiano-Obukhov scaling, respectively. Here, the nondimensional part is 
denoted by C^. The dimensional constants e, Se, a and g are the energy dissipation rate, the entropy dissipation 
rate, the thermal expansion coefficient and the gravitational acceleration, respectively. 

In order to determine the value of C^/A, we use the exit time statistics. The exit time, TE{r; p), is the time it takes 
for the separation of the particle pair to reach the threshold of pr from that of r. Assuming that the PDF of the pair 
separation obeys Eq. ([3|), the following expression of the mean exit time, 

A s(2s - 6) 

is derived^^ in the same manner as Ref. 14. The proportionality of the mean exit time to A^^i^p" — l)r* is also 
understood from the scaling law for the characteristic time, Eq. (|2ap . If we set the values of s, b and p and estimate 
that of the proportionality coefficient of in the right hand side of Eq. (|Bip from the data of DNSs or experiments, 
then the value of AjX can be calculated. Thus, we obtain the value of C^/A dividing AjX by the dimensional part of 
A, corresponding to the scaling of the DNS or experiment. 
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